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l/^ ■ ABSTRACT 

, The unexpected rising flux of early-type galaxies at decreasing ultraviolet (UV) wave- 

lengths is a long-standing mystery. One important observational const raint is the cor- 
relatio n between UV— optical colours and Mg2 line strengths found by iBurstein et al.l 
(|1988[) . The simplest interpretation of this phenomenon is that the UV strength is 
related to the Mg line strength. Under this assumption, we expect galaxies with larger 
• • I Mg gradients to have larger UV colour gradients. By combining UV imaging from 

. GALEX, optical imaging from MDM and SAURON integral-field spectroscopy, we in- 

k> ' vestigate the spatially-resolved relationships between UV colours and stellar popu- 

j_j I lation properties of 34 early-type galaxies from the SAURON survey sample. We find 

i that galaxies with old stellar populations show tight correlations between the FUV 

colours (FUV— and FUV— NUV) and the Mgb index, H/3 index and metallicity 
[Z/H]. The equivalent correlations for the Fe5015 index, a-enhancement [a/Fe] and 
age are present but weaker. We have also derived logarithmic internal radial colour, 
measured line strength and derived stellar population gradients for each galaxy and 
again found a strong dependence of the FUV— F and FUV— NUV colour gradients on 
both the Mg b line strength and the metallicity gradients for galaxies with old stellar 
populations. In particular, global gradients of Mg5 and [Z/H] with respect to the UV 
colour (e.g. A(Mg5)/A(FUV-NUV) and A[Z/H]/A(FUV-NUV)) across galaxies are 
consistent with their local gradients within galaxies, suggesting that the global cor- 
relations also hold locally. From a simple model based on multi-band colour fits of 
UV upturn and UV-weak galaxies, we have identified a plausible range of parameters 
that reproduces the observed radial colour profiles. In these models, the centers of 
elliptical galaxies, where the UV flux is strong, are enhanced in metals by roughly 
60% compared to UV-weak regions. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
photometry - galaxies: stellar content - galaxies: structure - ultraviolet: galaxies 
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1 INTRODUCTION 

Far ultraviolet (FUV) radiation was first detected in early- 
type galaxies (NGC 1291, 1316, 1553, 4406, 4486, 4649 
and the bulge of M31 ) by the Orbiting Astronomical 
Observatory-2 in 1969 (jcodel 1 19691 : ICode. Welch fc Pagel 
1 19721 '!. This discovery was surprising, since it had been ex- 
pected that such old populations would be entirely dark in 
the FUV. Even more surprisingly, a few elliptical galaxies 
showed extended rising flux for decreasing wavelengths be - 
low « 2000 A (e.g. lCode fc Welchll 19791 : iBertola et aLlll980D . 
This was called the UV upturn phenomenon. 

Early-on, one of the favored mechanisms for the UV up- 
turn was low-level star formation. Indeed, recent UV imag- 
ing surveys from the Galaxy Evolution Explorer (GALEX) 
have shown that a very large fraction (at least 15 %) of early- 
type galaxies as classified by the Sloan Digital Sky Survey 
(SPS S) exhibit sign a tures of ongoing or recent star forma- 
tion (|Yi et al.1 l2005l : ISchawinski et al.1 l2007l : | Kavirai et al.l 
[2003). In Ultraviolet Imagin g Telescope (UIT) images of 
the centre of M31, however. iBohlin et al.l |l98^ had not 
detected point sources suggesting that there is no evidence 
for the presence of main-sequen ce O stars ev en though its 
UV colour is blue. Furthermore, IWelchI 11982) claimed that 
the absence of strong CIV and Si IV in the FUV spectra of 
M31 indicates that young stellar populations are not present 
in the bulge of this galaxy. 

What, then, is the origin of the UV upturn phe- 
nomenon? Early stellar population models reproducing 
the red colours and spectral energy distributions (SEDs) 
of early-type galaxies in the optical band suggested 
that most stars in elliptic a l gal a xies should b e old 
and metal-rich (e.g., iFabeij 1 19721 : iTinslevI 1 19721 ). but 
these models failed to reproduce the strong observed 
UV flux. It is now generally recognized that early-type 
galaxies also contain old stellar components with tem- 
peratures exceeding 20,000 K. The leading hypothesis 
suggests that the FUV flux originates from a minor- 
ity population in the helium-burning horizontal-branch 
(HB) phase and its progeny, but there are two compet- 
ing H B solutions : met al- poor |Lee. Demaraue fc ZinnI 
19941 : IPark fc Leel 1 19971 : H uzzoni fc Gonzal cz-Lopczlira 
and m etal-rich jGreggio fc Rcnzini 



Horch. Demaraue fc Pinsonneaultl 



iBressan. Chiosi fc Fa.gottd 1 19941: 

Dorman. O'Connell fc Rood Il995l : lYi. Demaraue fc KimI 
19971 ). It is widely known that metal-poor HB stars can 
be hot and become UV bright when they are old. The 
oldest stars in this scenario are the most metal-poor and 
are most likely to be found in galactic centres, where the 



jYi et all 199£ 


). On 


(iBressan et al.l 


1994 



lYi. Demaraue fc Kim|[l997l ) generate UV sources by skip- 
ping the asymptotic giant branch (AGB) phase and require 
a similar age as the Milky Wa y halo. This phenomenon , 
dubbed the AGB manque stage (jGreggio fc Renzini|[l990h . 
is more pronounced at high values of helium abundance 
iHorch et all Il992l : [Porman et al] Il995h . Metal-rich stars 
thus become UV-bri ght sources . A tho rough review of these 
issues IS available in lO'Connell (|l999l ). 

Based on recent observations, a number of studies have 



identified peculiar globular clusters with extended hori- 
zontal branch (EHB) stars attributable to the presence 
of super-helium-rich populations (see e.g. I Lee et aD l2005l : 
iLee. Gim fc Casetti-Dinescull2007h . These studies also dis- 
cussed the discovery of numerous UV-bright gl obular clus- 
ters in the giant elliptical galaxy M87 (see e.g. ISohn et al.l 
l2006h . Although the origin of this helium enhancement is 
not yet fully understood, it could be a source of the FUV 
flux i n elliptical g alaxie s. Binary stars offer another possi- 
bility: iHan et al.1 (|2003f ). who constructed population syn- 
thesis models including binaries, concluded that most UV 
sources in elliptical galaxies come from binary stars and 
tha t there is no temporal evolution in FUV— V colour (see 



e^ 



Maxted et ahllioOll : IHan. Podsiadlowski fc Lvnas-Gravl 



l2007i r The other important candidate is post-AGB (PAGB 



stars (e.g. Greggio fc Renzinil Il990l : iDormari et al.1 [l9& 
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iBrown et al.lll997l . 120001 ). which also generate FUV radia- 
tion, but are thought to account for only 10 to 30 % of the 
UV flux. 

A crucial observational result is the dependence of the 
UV upturn on the stellar population and dynamical proper- 
ties of ea rly-type gal a xies jFabei 1 19831 : iBurstein et al.lll988l . 
hereafter IBBBFL88I ). lBBBFL8a found correlations of the 
UV— optical colours with both the Mg2 index and the cen- 
tral velocity dispersion, based on 24 quiescent early-type 
galaxies. However, a correlati on opposite to the Burstein re- 
lation has also been reported. iRich et al.l (|2005l ) constructed 
a FUV-r colour from GALEX UV and SDSS optical imag- 
ing and reported no correlation between the FUV-r colour 
and the Mg2 index, even after attempting to exclude galaxies 
with an active galactic nucleus (AGN) and star formation. 

Almost all previous studies focused on integrated UV- 
line strength relations. However, it is important to look at 
these relations as a function of radius within galaxies, since 
this can provide a finer discrimination between UV upturn 
theories. It has been known that a feature of the UV up- 
turn is the strongest in the centers of early-type galaxies. 
This produces the positive UV colour gradient that can be 
interpreted in ter ms of the difference in the stellar pop - 
ulations (see e.g. lO'ConneU et all 1 19921 : lO'Connelll 1 1 9991 ). 
iBrown et all (| 19971 )" claimed that the FUV light originates 
in a population with high metallicity and enhanced helium. 
If the UV colour variation is indeed related to the metal 
abundance, galaxies with larger metallicity gradie nts should 
have larger UV colour gradients. lOhl et aL 1 (119981 ) found no 
correlation between the FUV— -B colour gradients and inter- 
nal metallicity gradients via the Mg2 index, but they had 
only a small sample. 

To reconcile some of these discrepancies, we ob- 
tained spatially-resolved galaxy data that can provide de- 
tailed information on the stellar population distribution. 
By combining UV imaging observations from GALEX 
with SAURON integral-field early-type galaxy spectroscopy, 
we re-inve stigate the relationshi ps originally suggest ed by 
IBBBFL88I. In a com panion paper. [Bureau et all (|201ll . here- 
after IPaper XVIIll) used GALEX and SAURON data to re- 
visit the integrated UV-line strength relations using identi- 
cal apertures for all quantities. They recovered correlations 
of the integrated F\JV-V and FUV-NUV colours with the 
integrated Mg6 line strength index. Furthermore, they ar- 
gued that most outliers are due to galaxies exhibiting low- 
level star formation. 
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In this paper, we focus on the internal gradients of 
colour, measured line strength and derived stellar popula- 
tion property within each galaxy. In Section O we present 
a brief summary of the photometric and spectroscopic data 
and related analyses. The UV-index relations are presented 
and discussed in Sections [3] and 3] In Section [S] we discuss 
the origin of the UV upturn phenomenon based on several 
population synthetic models. We summarize our results and 
discuss their implications in Section [S] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 UV observations 

We observ ed 34 early-type galax ies from the SAURON sam- 
ple of 48 (|de Zeeuw et al.l |2002| ) with the medium-depth 
imaging mode of GALEX in both FUV (1350-1750 A) and 
NUV (1750-2750 A), as part of our own UV imaging sur- 
vey of the SAURON sample {GALEX guest investigator pro- 
grammes GIl-109 and GI3-041) and the GA LEX Nearby 
Galaxy Survey (NGS: iGil de Paz et al.ll2007l ). The typical 
exposure time per field was one orbit (« 1700 s), sufficient 
to exploit the two-dimensional nature of the images. Of 
the 48 early-type galaxies in the SAURON sample, another 
7 galaxies have only short exposures (typically 100-150 s) in 
the GALEX All-sky Imaging Survey (AIS), too shallow for 
spatially-resolved work. The remaining 7 galaxies have not 
been observed with GALEX, explaining the current sam- 
ple of 34 as listed in Table [2] Details of the data reduc- 
tion and analysis for these 34 early -type gal axies are de - 
scribed more fullv lJeong et al.l l|2009ll . hereafter Faper XII J) , 
so we only summarise our processing here. Details of the 
G ALEX instruments, pipeli ne and calibration are descri bed 
in lMartin erall (|2005l ') and lMorrissev et all (|2005l. 120071 ') . 

The spatial resolution of the GALEX images is approxi- 
mately 4'.'5 and 6'.'0 FWHM in FUV and NUV, respectively, 
sampled with l'.'5 x l'.'5 pixels. Although the images were 
delivered after pre-processing, we performed our own sky 
subtraction by measuring the sky level in source-free re- 
gions of the images. Additionally, we convolved the FUV 
data to the spatial resolution of the NUV observations be- 
fore any analysis, to avoid spurious colour gradients in the 
central regions. We also generated a mask image to min- 
imize contamination by nearby sources using SEXTRAC- 
TOR, and carefully examined. The mask images are good 
in most cases, but for a few cases, we had to mask some 
sources manually. In the case of NGC 4486, for example, we 
mask the northwestern quadrant to minimize the effect of 
the jet, but we do not remove the likely contamination by 
non-thermal UV emission. Finally, we carried out surface 
photometry by measuring the surface brightness along el- 
liptical annuli using the ELLIPSE task within the STSDAS 
ISOPHOTE package in IRAF (Image Reduction and Anal- 
ysis Facility), after masking out nearby sources. For each 
image, the centre of the isophotes was fixed to the centre of 
the light distribution and the position angle (PA), elliptic- 
ity (e) and surface brightness (fj,) were fitted as a function 
of the radius. The ellipses were fitted only to the NUV im- 
ages, which have much better signal-to-noise ratios (S/N) 
at all radii than the FUV images, and these ellipses were 
then superimposed on the FUV images to derive meaning- 
ful colours. We also de-reddened the colours for Galactic 



extinction using values of ^fuv = 8.376 x E(B— V) and 
^Nuv = 8.741 X E(B- V) jWvder et aLll2005l ) and the red- 
dening maps of lSchlegel. Finkbeiner fc David (|l998l ). 



2.2 Optical observations 

Ground-based optical imaging observations with the Hubble 
Space Telescope (HST) filter F555W (similar to Johnson V^) 
were obtained using the MDM Observatory 1.3-m McGraw- 
Hill Telescope over several observing runs in 2002-2005, 
again as part of a larger survey targeting the whole SAURON 
galaxy sa mple. Details of the MDM observat ions are de - 
scribed in lFalcon-Barroso et al.l (|201ll . hereafter lPaper XIXl) 
and the data were reduced and calibrated in the standard 
manner. The field-of-view of the MDM images is 17^3 x 17^3 
with 0'.'508 X 0'.'508 pixels, allowing for accurate sky subtrac- 
tion and proper sampling of the seeing. The seeing for the 
observations was ~ l'.'2. 

The MDM data were also convolved to the resolution of 
the GALEX NUV data for the current analysis. The surface 
brightness profiles were measured after superimposing the 
NUV ellipses as for the FUV data. The magnitudes at each 
isophote were then corrected for Galactic extinction using 
a value of Rv = 3.1 (ICardelh. Clavtori fc Mathislll989h and 
the reddening maps of Schlegel et al. ( 19981 ). 



2.3 SAURON linestrengths and kinematics 

The SAURON observations were designed to determine 
the two-dimensi onal stella r kinematics ([Emsellem et al.l 
|2004 hereafter IPaper Ilj: lEmsellem et al.l |2007| . here- 
after iPaper l5([). stellar line strengths a nd popu- 
lations (iKuntschner et al.l |2006|. hereafter IPaper vj 
iKuntschner et aU 2010l . hereafter Paper XVIll) and ionised 
gas kinematics and properties of 48 nearby early- typ e 
galaxies in the field and clusters (|de Zeeuw et al.l 20021 ) . 
using the panoramic integral-field spectrograph for the 
William Herschel Telescope (WHT) based on the SAURON 
TIGER microlens concept . The observations and data 
reduction are described in the indicated papers. 

To explore the UV-line strength relations as a function 
of radius within individual galaxies, we also imposed the el- 
lipses from the NUV images on the SAURON maps of measured 
line strengths (Mg&, Fe5015 and II/3), derived (as opposed 
to measured) stellar population properties (a-enhancement 
[a/Fe], metallicity [Z/ H] and age) and stellar velocity disper- 
sion (a). iPaper XVlJ transformed the Mgfo, Fe5015, Fe5270 
and II/3 line strengths into stellar population prop erties us- 
ing th e single stellar population (SSP) models of SchiavonI 
(|2007l ). This process is fully described in iPaper XVlj . It"7s 
important to bear in mind that these parameters are de- 
rived from only three spectral lines (sometimes four lines), 
owing to the narrow bandwidth of the SAURON spectra, and 
that they are "SSP-equivalent" parameters, relying on the 
supposition that the stars in each galaxy have formed in a 
single burst. This is an ad hoc assumption; however, it al- 
lows comparisons between the SSP-equivalent values of dif- 
ferent sample galaxies. We use simple luminosity- weighted 
line strength, population and velocity dispersion averages for 
each isophote . The legitimacy of this method is discussed in 
IPaper XVIlS 
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Figure 1. UV upturn dependence on measured line strengths and derived stellar population properties. The Mgfe {top-left), Fe5015 
{middle-left), H/3 {bottom-left), [a/Fe] {top-right), [Z/H] {middle-right) and age {bottom-right) parameters arc shown as a function of 
the FUV— y colour. We note that all line strengths and derived stellar population properties are based on the SAURDN data published in 
Paper VI and Paper XVII, respectively. The sample is divided into quiescent early-type galaxies (circles, 11/3^^/2 ^ 1-8 ^) s-i^d galaxies 



with recent star formation (stars, 



:o/2 



> 1.1 



The solid line in each panel is a linear fit to all of the quiescent galaxy profiles. 



weighting each radius by its V-band luminosity, and the correlation coefficient is reported in the bottom-right corner. The symbols are 
colour-coded according to the integrated stellar velocity dispersion within one effective radius (ce ) . All the points for an individual galaxy 
are represented by a single colour, and the symbol for the central point is larger. 
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Table 1. Parameters of the best-fit Hnear UV-linc strength (and stellar property) relations for old galaxies only. 



v_yL/iW HI 


Index 


Slope 




Unit 




Unit 






FVY-V 


Mgfe 


-0.48 ± 0.03 


A 


mag~^ 


7.97 ±0.23 


A 


-0.77 


<0.001 




Fe 


-0.18 ±0.04 


A 


mag~^ 


6.33 ±0.29 


A 


-0.43 


<0.001 




H/3 


0.13 ±0.01 


A 


mag~^ 


0.58 ±0.09 


A 


0.76 


<0.001 




[a/Fe] 


—0.04 ± 0.01 




mag 


0.55 ± 0.05 




-0.46 


<0.001 




[Z/H] 


-0.12 ±0.02 




mag~^ 


0.79 ± 0.08 




-0.59 


<0.001 




age 


-0.06 ±0.01 


Gyr 


mag^^ 


1.56 ±0.05 


Gyr 


-0.49 


<0.001 






-0.14 ±0.02 


km s"-"^ 


mag^^ 


3.15 ±0.09 


km s^^ 


-0.54 


<0.001 


FUV-NUV 


Mgb 


-0.73 ± 0.04 


A 


mag~^ 


5.06 ± 0.05 


A 


-0.86 


<0.001 




Fe 


-0.36 ± 0.06 


A 


mag~^ 


5.09 ±0.07 


A 


-0.59 


<0.001 




H/3 


0.19 ±0.02 


A 


mag^-'^ 


1.39 ±0.03 


A 


0.72 


<0.001 




[a/Fe] 


-0.03 ±0.01 




mag^^ 


0.33 ±0.01 




-0.29 


0.005 




[Z/H] 


-0.20 ±0.02 




mag^-'^ 


0.18 ±0.02 




-0.72 


<0.001 




age 


-0.08 ±0.01 


Gyr 


mag~^ 


1.21 ±0.01 


Gyr 


-0.50 


<0.001 






-0.25 ±0.03 


km s^-"^ 


mag^^ 


2.50 ±0.02 


km s^^ 


-0.68 


<0.001 



3 UV RELATIONS 

As argued in the introduction, spatially-resolved galaxy data 
provide valuable information for discriminating between UV 
upturn theories. The internal stellar population gradients 
within each galaxy can yield their own "local" relations, 
rather than relying on UV-stellar population correlations 
from the inte grated measurements of many galaxies (see e.g. 
IPaper XVIll . When exploring correlations between stel- 
lar population parameters and multi-band photometry, our 
main advantage is that all photometric and spectroscopic 
measurements are derived using the same ellipses. We note 
that all stellar absorption line strengths (Mg6, Fe5015 and 
H/3), stellar population properties {[a/Fe], [Z/H] and age) 
and stell ar kinematics (a) are based on the SAUR ON data pu b- 
Ushed in IPaper Vj. IPaper XVlj IPaper Hj and IPaper IXl . 

Strong Balmer absorption lines betray the pres- 
ence of young stars. Based on the prediction of 
various stellar population synthesis models ( see e.g. 
iThomas. Maraston fc Bended [20o3 : ISchiavonI 120071 ). an H/3 
line strength of « 1.8 A can be used to divide the sam- 
ple into galaxies with effectively only old stellar popula- 
tions and those which have a contribution from younger 
stars. To minimise the effects of galaxies with recent star 
formation on the correlations, we represent galaxies with 
an integrated H/3 line strength withi n half an effec tive ra- 
dius H/3jj^/2 > 1.8 A (see Table 1 in iPaper XVllt as star 
symbols (Figures 1-7, 10 and 11). The symbols in Figures 
1-7, 10-12 are also colour-coded according to the integrated 
stellar veloc ity dispersio n within one effective radius (ae; 
Table Al in lPaper XVIll) . descending from red to purple. 



3.1 Radial profiles 

The left-hand panels of Figure [T] show the correlations be- 
tween the FUV— colour and the line strength indices Mg6, 
Fe5015 and H/3 as a function of radius. Linear fits to all the 
quiescent galaxy (H/3jj^/2 ^ 1-8 A) profiles, after weighting 
each radius by its luminosity at V band, are shown as solid 
lines. The value of the correlation coefficient is reported in 
the bottom-right corner of each panel, and the parameters 
of the best-fit lines are reported in Table[T] with errors prop- 
erly taken into account. 



The top-left panel of Figure [T] shows the (FUV-V)- 
Mgb "Burstein" relation. This relation shows the tightest 
correlation among the different line strength indices, as dis- 
cussed in many previous studies. The slope of a linear fit to 
the integrated properties of the quiescent galax ies (including 
the 7 AIS galaxies; see Table 1 of lPaper XVH^ is -0.40±0.09 
A mag"'^, whereas the slope measured here based on the 
luminosity-weighted galaxy profiles (solid line; global gradi- 
ent) has a value of -0.48±0.03 A mag . These two slopes are 
thus consistent within th e errors and o ur results are in good 
agreement with those of lPaper XVIlj . All of the quiescent 
galaxies exhibit bluer FUV— 1^ colours at smaller radii, sug- 
gesting that the FUV light is more centrally-concentrated 
than the optical fight, also in agreement w ith earlier studies 
(see e.g. lOhl et al.|[l99l : IO'Connell|[l999l '). 

We stress that galaxies with young stellar populations 
(H/3fl^/2 > L8 A, star symbols) show different trends (either 
opposite trends or constant Mg& values as FUV— V colours 
are varied) than quiescent galaxies, so the relation becomes 
less distinct. These galaxies also have lower Mgb values. It is 
clear therefore that a tight (FUV— V^)-Mg6 relation is only 
present for galaxies which show no sign of on-going or recent 
star formation. 

The middle-left panel of Figure [T] shows the relation 
of FUV— y with Fe 5 015. T his r elation is mu ch weaker, as 
found by lO'Connelll j 19991 ) and IPaper XVIll This implies 
that the UV strength is more strongly linked to the behavior 
of the lighter elements (e.g. Mg), not the iron peak. Q We 
also note that galaxies with a young population appear at 
lower Fe values. 

In the case of H/3 (bottom-left panel of Figure[T]), there 
is a clear correlation with the FUV— 1/ colour. A tight corre- 
lation is, again, only present among quiescent galaxies. The 
separation of young galaxies (stars) from old ones (circles) is 
however somewhat artificial in this plot, as we use H/3 itself 
as the dividing criterion. The observed variation of H/3 with 
FUY—V is also expected from the (FUV— V')-Mg b corre- 
lation. According to stellar population models, differences 

^ Admittedly, our analysis is limited to Fe5015 and Mgb as iron- 
peak and magnesium strength indieators due to the narrow wave- 
length coverage of the SAURON instrument. A wider wavelength 
coverage is desired to constrain them more robustly. 
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Figure 2. Same as Figure [T] but for tlie NUV— V colour. 



in metallicity cause variati ons in the H/3 l i ne str ength. For 
example, in the models of I Thomas et al.l (|2003l ) , an Mg b 
value of 5.16 A and an age of 12 Gyr correspond to an 
H/3 line strength of 1.50 A, whereas for Mgfe = 3.11 A and 
age = 12 Gyr the H/3 line strength increases to 1.77 A. The 
observed trend can thus easily be explained by the effect of 
metallicity on H/3. 

The right-hand panels of Figure [T] show the correlations 
between the FUV— F colour and the stellar population prop- 



erties ([a/Fe], [Z/H] and age) as a function of radius. The 
parameters of the best-fit lines and the values of the cor- 
relation coefficients are again listed in Table [1] Only [Z/H] 
shows a modest correlation with FUV. It is surprising that 
[a/Fe] does not show a clear correlation, considering the 
tight correlation of Mg6 with FUV. 

Analogously to Figure[T] Figure[2]shows the correlations 
between the NUV— V colour and the line strength indices 
(Mg6, Fe5015 and H/3) and stellar population parameters 
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Figure 3. Same as Figurc[T]but for the FUV-NUV colour. 



([a/Fe], [Z/H] and age) as a function of radius. In contrast 
to Figure [TJ we note no correlation between the NUV— V 
colour and those quantities. The UV upturn phenomenon 
clearly dominates the spectral energy distribution only at 
wavelengths blueward of 2000 A. We therefore believe that 
these trends are caused by two main effects: young stars 
and line blanketing. Recent star formation {t < 1 Gyr) has 
a greater effect in the NUV than in the FUV, as only on- 
going star formation {t < 0.1 Gyr) can contribute to the 



FUV light. On the other hand, line blanketing is signifi- 
cant in the NUV as metallicity increases, as a large number 
of partially overlapping absorption lines can suppress the 
continuum. These are probably the main reasons why the 
FUV correlations become weaker in the NUV. Nevertheless, 
galaxies with extremely high Mg b line strength and a strong 
UV upturn also show bluer NUV— V" colours. 

Figure[3]shows the correlations between the FUV-NUV 
colour and the line strength indices (Mgb, Fe5015 and H/3) 
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Figure 4. Same as Figurc[T]but for the FUV— NUV and FW—V colours as a function of the integrated stellar velocity dispersion within 
one effective radius (To. 



and stellar population parameters ([a/Fe], [Z/H] and age) 
as a function of radius. The luminosity-weighted linear fits 
to the data with H/3fl^/2 ^ 1.8 A are shown as solid 
lines and the parameters of the best-fit lines are again re- 
ported in Table [T] The value of the correlation coefficient is 
also reported in the bottom-right corner of each panel and 
Table [T] Note that the UV-line strength correlations are 
much more pronounced when the FU V— V colour is repl aced 
by the FUV- NUV colour (see also IPonas et al.ll2007l and 
IPaper XVIl3) . For example, the (FUV-NUV)-Mgfe corre- 
lation is tighter than that of (FUV-y)-Mgb. We men- 
tioned above that line blanketing is substantial in the NUV 
and increases with metallicity. This is likely the reason why 
FUV-NUV correlations are much tighter than those with 
FUV— Another advantage of using FUV-NUV rather 
than FUV-V is that both FUV and NUV magnitudes are 
measured with the same instrument, and uncertainties in 
the UV to optical flux calibration are avoided. 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 
A(FUV-NUV)/Alog(R) (This paper) 



3.2 Stellar velocity dispersion 

It is well-known that there is a tight correlation between 



e.E. iTerlevich et al.l 


19811: iBender. Burstein & Fabeij Il993l: 


ICoUess et al.l 19991: 


Bernardi et al. 20031). This is a well- 



studied relation linking a dynamical quantity a with a quan- 
tity (the Mg index) that depends on the stellar populations. 
A number of studies have implied a correlation of the UV- 
upt urn phenomenon with the stellar ve locity dispersion (see 
e.g. lBBBFL8lll988l : lDonas et al.ll2007l ). a good measure of 
a galaxy's gravitational potential depth and total dynamical 
mass. 

Figure [4] shows the colour-tJo relations. Comparing to 
Figures [T] and [3] (top-left), one can see that the (FUV— F)- 
Mg6 and (FUV-NUV)-Mg& relations are tighter than the 
(FUV-F)-o-c and_(FUV-NUV)-cre relations, confirming 
the results of lBBBFL8l . 



Figure 5. Comp arison of our (FUV- NUV) radial colour gradi- 
ents with those of lCarter et al. 1 1I2OIII) . 



3.3 Internal gradients 

From the radial surface brightness, line strength and stel- 
lar population profiles, we derive logarithmic internal ra- 
dial gradients for each galaxy using least-square fits. The 
radial range of the fits generally extends from the FWHM 
of the UV point spread function to half the effective ra- 
dius (i?e/2). The radial gradients are thus defined by their 
global slopes. In addition, to estimate the gradient uncer- 
tainties, we performed Monte-Carlo simulations comprising 
10'^ realisations, adding noise compatible with the observa- 
tions to each data point. In Figure [S] we show the compar- 
is on of our (FUV-NU V) radial colour gradients with those 
of lCarter et al.l (|201ll ) . The agreement is generally good: the 
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Table 2. Logarithmic internal radial gradients. 



/~,„i„„„ AFUV-V AFUV-NUV AMgb AFc5015 AHb A[a/Fc] A[Z/H] Alog(agc) 

vjdiaxy Alog(R) Alog(R) Alog(R) Alog(R) Alog(R) Alog(R) Alog(R) Alog(R) 

(mag dex^-"-) (mag dex^-"-) (A dex"-"-) (A dex^-"-) (A dcx"-"-) (Gyr dex^-^) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 
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Columns: (1) Galaxy identifier. (2)-(3): Radial colour gradients. (4)-(9): Radial line strength and stellar population gradients. 



mean of the absolute value of the (FUV— NUV) radial colour 
gradient difference is ~0.55 magdex"^. The largest differ- 
ence is found for NGC 4552 likely due to the determination 
of the radial range of the fits. We adopted the same radial 
range to derive the internal radial gradients for each galaxy. 
The logarithmic internal radial gradients of the UV colours, 
line strengths and stellar population properties are listed in 
Tabled 

We plot the internal gradients of the line strength in- 
dices and the stellar population parameters as a function of 
the FUV— 1/ internal colour gradients in Figure |6l To pre- 
vent any dilution of the correlation attributable to galaxies 
with recent star formation, we exclude galaxies with evi- 
dence of young stars (H/3/i^/2 > 1-8 A) from the fit and plot 
them as grey stars. The value of the correlation coefficient is 
reported in the bottom-right corner of each panel, and the 
parameters of the best-fit lines are reported in Table [S] 

Most importantly, only Mg6 and [Z/H] internal gradi- 
ents for quiescent galaxies show correlations. This suggests 
that either Mg is a primary driver of the UV upturn, given 
that Mg fe is a reliable tracer of Mg, or Mg and the UV up- 
turn share the same primary driver. If the stellar population 



properties ([a/Fe], [Z/H] and age) are reliable, we conclude 
further that the UV strength is not predominantly governed 
by the a-enhancement. The correlation between the [Z/H] 
index and the UV colour can be attributed to the role of Mg 
in the [Z/H] measurement. The weak or absent correlation 
with a-enhancement ([a/Fe]), however, is contrary to our 
expectations. 

Analogously to Figure [S] Figure [7] shows the 
FUV— NUV logarithmic internal radial colour gradients ver- 
sus the internal gradients of the line indices and stellar 
population properties. The (FUV— NUV)-Mg6 correlation 
is again much tighter than the (FUV— y)-Mg6 correlation, 
as discussed in Section [3TT] and the correlation with Mgb is 
much stronger than that with [Z/H]. 

Figure |8] shows the internal gradients of the velocity 
dispersion as a function of the UV internal colour gradients. 
Comparing to Figures |5] and [3 the correlations with are 
much weaker than those with both Mgfo and [Z/H]. This 
suggests that the UV upturn is driven by stellar population 
properties rather than the velocity dispersion. 

For clarity, we divide galaxies into groups exhibiting 
similar UV radial colour profiles for the region interior to 
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Table 3. Parameters of the best-fit Unear radial gradients for old galaxies only. 



Colour 


Index 


Slope 


Zero- point 


Coefficient 


Significance level 


FUV-y/Alog(R) 


Mgfe/Alog{R) 


-0.35 ± 0.13 


-0.22 ± 0.12 


-0.59 


0.001 




Fe/Alog(R) 


-0.21 ± 0.16 


-0.30 ±0.15 


-0.31 


0.14 




H/3/Alog{R) 


-0.01 ± 0.09 


0.01 ±0.07 


-0.07 


0.67 




[a/Fe]/Alog(R) 


— O.Ul ± 0.U5 


0.01 ± 0.05 


-0.07 


0.68 




[Z/H] /Alog{R) 


-0.08 ± 0.04 


-0.09 ± 0.04 


0.61 


<0.001 




log age/ Alog(R) 


0.02 ± 0.04 


-0.01 ± 0.04 


0.21 


0.28 




log<7e/Alog{R) 


-0.30 ±0.04 


-0.12 ± 0.04 


-0.21 


0.27 


FUV-NUV/Alog{R) 


Mgfe/Alog{R) 


-0.71 ± 0.16 


0.10 ±0.14 


-0.83 


<0.001 




Fe/Alog(R) 


-0.31 ± 0.25 


-0.20 ±0.23 


-0.31 


0.14 




H/3/Alog{R) 


0.13 ± 0.13 


-0.10 ±0.11 


0.44 


0.018 




[a/Fe]/Alog(R) 


-0.07 ±0.08 


0.05 ± 0.07 


-0.28 


0.18 




[Z/H] /Alog{R) 


-0.11 ± 0.06 


-0.06 ±0.06 


0.54 


0.004 




log age/ Alog(R) 


-0.03 ± 0.06 


0.04 ±0.06 


-0.20 


0.31 




log(Tc/Alog{R) 


-0.08 ± 0.06 


-0.09 ± 0.05 


-0.30 


0.15 



Table 4. Internal gradients. 



p , AlogcT Alogg AMgb AMgb A[Z/H] A[Z/H] UVt™ 

^ AFUV-V AFUV-NUV AFUV-V AFUV-NUV AFUV-V AFUV-NUV ^ ^ >^ 

(dex mag~^) (dex mag"-"^) (A mag"-"^) (A mag"-'^) (dex mag"-"^) (dex mag"-') 

(1) (2) (3) (4) (5) (6) (7) (8) 



NGC0474 


0, 


,28 


± 


0, 


,15 


-0, 


,22 


± 


0.07 


2, 


,59 


± 


0, 


,72 


-1, 


,94 


± 


0, 


,32 


0, 


,26 


± 


0, 


,30 


-0, 


,15 


± 


0, 


,14 


RSF 


NGC0524 


-0, 


,14 


± 


0, 


.05 


-0, 


.19 


± 


0.05 


-0, 


.40 


± 


0. 


,24 


-0, 


.47 


± 


0, 


,25 


-0, 


,02 


± 


0, 


,10 


0, 


,02 


± 


0, 


.10 


UV-weak 


NGC0821 


-0, 


,09 


± 


0, 


,07 


-0, 


,05 


± 


0.04 


-0, 


.49 


± 


0, 


,31 


-0, 


,28 


± 


0, 


,17 


-0, 


,17 


± 


0, 


,14 


-0, 


,09 


± 


0, 


,08 


UV-weak 


NGC1023 


-0, 


,79 


± 


0, 


,16 


-0, 


,50 


± 


0.07 


-0, 


.48 


± 


0, 


,74 


-0, 


,18 


± 


0, 


,31 


-0, 


,14 


± 


0, 


,31 


-0, 


,11 


± 


0, 


,13 


RSF 


NGC2695 


-0, 


,28 


± 


0, 


.05 


-0, 


.26 


± 


0.04 


-0, 


.20 


± 


0. 


,22 


-0, 


.20 


± 


0, 


,21 


-0, 


.04 


± 


0, 


.09 


-0, 


,04 


± 


0, 


.08 


UVX 


NGC2699 


-0. 


,24 


± 


0, 


.04 


-0, 


.18 


± 


0.03 


-0, 


.64 


± 


0. 


,18 


-0, 


.49 


± 


0, 


,14 


-0, 


.47 


± 


0, 


.08 


-0, 


,35 


± 


0, 


.06 


UV-weak 


NGC2768 


0. 


,04 


± 


0, 


.32 


-0, 


.13 


± 


0.14 


2, 


.39 


± 


1, 


,52 


-1, 


.50 


± 


0, 


,66 


1, 


,21 


± 


0, 


.66 


-0, 


,54 


± 


0, 


.29 


UV-weak 


NGC2974 


-0, 


,12 


± 


0, 


,04 


-0, 


,13 


± 


0.04 


-0, 


,46 


± 


0, 


,17 


-0, 


,47 


± 


0, 


,18 


-0, 


,14 


± 


0, 


,07 


-0, 


,14 


± 


0, 


,07 


RSF 


NGC3032 


-0, 


,04 


± 


0, 


,01 


-0, 


,09 


± 


0.02 


0, 


,49 


± 


0, 


,04 


1, 


,24 


± 


0, 


,10 


0, 


,26 


± 


0, 


,02 


0, 


,60 


± 


0, 


,04 


RSF 


NGC4150 


-0, 


,10 


± 


0, 


.02 


-0, 


.12 


± 


0.02 


0, 


.01 


± 


0. 


,08 


0, 


,01 


± 


0, 


,09 


0, 


,00 


± 


0, 


.03 


0, 


,00 


± 


0, 


.04 


RSF 


NGC4278 


-0, 


,10 


± 


0, 


.04 


-0, 


.09 


± 


0.04 


-0, 


.87 


± 


0. 


,18 


-0, 


.89 


± 


0, 


,17 


-0, 


.20 


± 


0, 


.07 


-0, 


,21 


± 


0, 


.07 


UVX 


NGC4374 


-0, 


,08 


± 


0, 


,05 


-0, 


,08 


± 


0.04 


-0, 


.83 


± 


0, 


,24 


-0, 


,75 


± 


0, 


,21 


-0, 


,22 


± 


0, 


,10 


-0, 


,20 


± 


0, 


,09 


UVX 


NGC4387 


-0, 


,12 


± 


0, 


,07 


-0, 


.14 


± 


0.04 


0, 


.40 


± 


0. 


,30 


0, 


,32 


± 


0, 


,19 


0, 


,05 


± 


0, 


.12 


0, 


,00 


± 


0, 


.08 


UV-weak 


NGC4458 


-0, 


,26 


± 


0, 


.21 


-0, 


.27 


± 


0.11 


-0, 


.04 


± 


0. 


,99 


0, 


,05 


± 


0, 


,53 


0. 


,58 


± 


0, 


.44 


-0, 


,22 


± 


0, 


.27 


UV-weak 


NGC4459 


-0, 


,06 


± 


0, 


,01 


-0, 


,12 


± 


0.02 


-0, 


,17 


± 


0, 


,05 


-0, 


,38 


± 


0, 


,10 


-0, 


.15 


± 


0, 


,02 


-0, 


,30 


± 


0, 


,04 


RSF 


NGC4473 


0. 


,01 


± 


0, 


,15 


0, 


,03 


± 


0.05 


-1, 


,33 


± 


0, 


,68 


-0, 


,84 


± 


0, 


,24 


-0, 


,54 


± 


0, 


,31 


-0, 


,32 


± 


0, 


,11 


UV-weak 


NGC4477 


-0, 


,75 


± 


0, 


.41 


-0, 


.55 


± 


0.20 


-0, 


.53 


± 


1, 


,92 


-0, 


.58 


± 


0, 


,96 


-0, 


.44 


± 


0, 


.82 


-0, 


,36 


± 


0, 


,41 


UV-weak 


NGC4486 


-0, 


,11 


± 


0, 


.04 


-0, 


.18 


± 


0.07 


-0, 


.60 


± 


0, 


,17 


-1, 


.08 


± 


0, 


,31 


-0, 


.19 


± 


0, 


.07 


-0, 


,33 


± 


0, 


,12 


UVX 


NGC4526 


-0, 


,08 


± 


0, 


.03 


-0, 


.17 


± 


0.10 


0, 


.08 


± 


0, 


,16 


-0, 


.53 


± 


0, 


,45 


-0, 


.15 


± 


0, 


.07 


0, 


,35 


± 


0, 


.18 


RSF 


NGC4546 


-1, 


,03 


± 


0, 


,23 


-0, 


,33 


± 


0.07 


-3, 


,03 


± 


1, 


,10 


-0, 


.93 


± 


0, 


,33 


-1, 


,03 


± 


0, 


,52 


-0, 


,33 


± 


0, 


,16 


UV-weak 


NGC4550 


0. 


,44 


± 


0, 


,11 


0, 


,26 


± 


0.06 


-0, 


,66 


± 


0, 


,48 


-0, 


,43 


± 


0, 


,27 


-0, 


,64 


± 


0, 


,22 


-0, 


,37 


± 


0, 


,13 


RSF 


NGC4552 


-0, 


,06 


± 


0, 


.03 


-0, 


.09 


± 


0.04 


-0, 


.34 


± 


0, 


,12 


-0, 


.49 


± 


0, 


,17 


-0, 


.09 


± 


0, 


.05 


-0, 


,13 


± 


0, 


.07 


UVX 


NGC4564 


-0, 


,07 


± 


0, 


.03 


-0, 


.08 


± 


0.04 


-0, 


.30 


± 


0, 


,14 


-0, 


.35 


± 


0, 


,17 


-0, 


.09 


± 


0, 


.06 


-0, 


,10 


± 


0, 


.08 


UVX 


NGC4570 


-0, 


,31 


± 


0, 


,07 


-0, 


,23 


± 


0.05 


-0, 


,60 


± 


0, 


,31 


-0, 


,40 


± 


0, 


,24 


-0, 


.14 


± 


0, 


.15 


-0, 


,11 


± 


0, 


.12 


UV-weak 


NGC4621 


-0, 


,08 


± 


0, 


.03 


-0, 


.09 


± 


0.03 


-0, 


.48 


± 


0, 


,12 


-0, 


.49 


± 


0, 


,13 


-0, 


.16 


± 


0, 


.06 


-0, 


,17 


± 


0, 


.06 


UVX 


NGC5198 


-0, 


,15 


± 


0, 


.04 


-0, 


.15 


± 


0.04 


-0, 


.39 


± 


0, 


,19 


-0, 


.39 


± 


0, 


,21 


-0, 


.13 


± 


0, 


.09 


-0, 


,14 


± 


0, 


.09 


UVX 


NGC5308 


-0, 


,41 


± 


0, 


,09 


-0, 


,33 


± 


0.07 


-0, 


,76 


± 


0, 


,44 


-0, 


,60 


± 


0, 


,34 


-0, 


,31 


± 


0, 


,20 


-0, 


,24 


± 


0, 


,15 


UV-weak 


NGC5813 


0. 


,04 


± 


0, 


,15 


0, 


,07 


± 


0.08 


0, 


,10 


± 


0, 


,70 


-0, 


,10 


± 


0, 


,36 


0, 


,59 


± 


0, 


,30 


0, 


,23 


± 


0, 


,15 


A/Z 


NGC5831 


-0, 


,12 


± 


0, 


,06 


-0, 


,08 


± 


0.04 


-1, 


,33 


± 


0, 


,30 


-1, 


,20 


± 


0, 


,20 


-0, 


,20 


± 


0, 


,13 


-0, 


,19 


± 


0, 


,09 


A/Z 


NGC5838 


-0, 


,27 


± 


0, 


.04 


-0, 


.28 


± 


0.04 


-0, 


.50 


± 


0, 


,18 


-0, 


.52 


± 


0, 


,18 


-0, 


.14 


± 


0, 


.07 


-0, 


,15 


± 


0, 


.08 


UV-weak 


NGC5845 


-0, 


,17 


± 


0, 


.05 


-0, 


.11 


± 


0.03 


-1, 


.44 


± 


0, 


,23 


-0, 


.96 


± 


0, 


,16 


-0, 


.30 


± 


0, 


.10 


-0, 


,20 


± 


0, 


.07 


A/Z 


NGC5846 


-0, 


,03 


± 


0, 


.05 


-0, 


.04 


± 


0.06 


-0, 


.54 


± 


0. 


,22 


-0, 


.65 


± 


0, 


,28 


-0, 


.19 


± 


0, 


.09 


-0, 


,22 


± 


0, 


.11 


UVX 


NGC5982 


-0, 


,14 


± 


0, 


,04 


-0, 


,13 


± 


0.04 


-0, 


,59 


± 


0, 


,20 


-0, 


,59 


± 


0, 


,20 


-0, 


.20 


± 


0, 


,08 


-0, 


,21 


± 


0, 


.08 


UVX 


NGC7457 


-0, 


,02 


± 


0, 


,03 


-0, 


,02 


± 


0.03 


0, 


,02 


± 


0, 


,14 


0, 


,02 


± 


0, 


,12 


0, 


,02 


± 


0, 


,05 


0, 


,02 


± 


0, 


,04 


RSF 



Columns: (1) Galaxy identifier. (2)-(7): Gradients of the UV eolour-index relations. (8): UV-optical radial colour profiles classification. 

RSF: recent star formation. UVX: UV upturn. 
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A(FUV-V)/Alog(R) (mag dex"') 



1 2 

A(FUV-V)/Alog(R) (mag dex"') 



Figure 6. Same as Figure [T] but for the FUV— V logarithmic internal radial eolour gradient versus the measured and derived stellar 
population gradients of the sample galaxies. The gradients of Mgb (top-left), Fe5015 [middle-left), H/3 (bottom-left), [o/Fe] (top-right), 
[Z/H] (middle-right) and age (bottom-right) are shown as a function of the FUV— V colour gradient. Quiescent galaxies (H/9^^/2 1-8 A) 
are shown as solid circles, while galaxies with a young population are shown as grey stars. The solid line in each panel is a fit to the 
quiescent galaxies, and the correlation coefficient is reported in the bottom-right corner. 



the effective radius (see Figure 1 of IPaper XIll ) and listed 



in Table |4l This is unlike the classification of IPaper XIIll 
that considered much larger radial ranges, because SAURON 
covers each galaxy only up to about one effectiv e radius. 
The classification differences between IPaper Xllj and this 
paper thus arise from different radial ranges considered. To 
select purely old galaxies, we first exclude g alaxies with 
recent star formation as listed in Table 1 of iPaper Xllj . 



Galaxies included in this category are NGC 474, 1023, 2974, 
3032, 4150, 4459, 4526, 4550 and 7457. We then consider 
the NUV — V" radial c olour profiles. As discussed in Sec- 
tion 3.3 of IPaper Xllj. NUV— y colours usually show neg- 
ative slopes due to the gradients of the underlying stel- 
lar populations, mainly in age and metallicity. This age- 
metallicity degeneracy is well-known and it is difficult to 
distinguish the effects of a small change in age from those of 
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Figure 7. Same as Figure [6] but for the FUV— NUV colour. 



a small change in metallicity. Therefore, we classify a galaxy 
as an age-metallicity degeneracy galaxy if its NUV— V pro- 
file is below 5.0 mag arcsec"^ around the effective radius 
only. These galaxies include NGC 58f3, 583f and 5845. 
Secondly, we consider the FUV— NUV and FUV— F pro- 
files inside the effective radius. If a galaxy has a central re- 
gion with FUV-NUV< 0.9 mag arcsec"^ and FUV-V < 
6.21 mag arcsec~^ then, following the new classification 
scheme of lYi et al.l |2oi3), we classify it as a UV upturn 
galaxy. These galaxies are NGC 2695, 4278, 4374, 4486, 
4552, 4564, 4621, 5198, 5846 and 5982. Finally, we classify 
the remaining galaxies as UV-weak galaxies. Galaxies in- 



cluded in this category are NGC 524, 821, 2699, 2768, 4387, 
4458, 4473, 4477, 4546, 4570, 5308 and 5838. 

Furthermore, we calculate the gradients of the UV 
colour-index relations (e.g. A(Mg6)/A(FUV— F)) for each 
individual galaxy using least-square fits (see Table [Jjl in the 
same way as described in Section [3.31 We note that the val- 
ues of e.g. [A(Mgb)/A(FUV-V)] are not equal to those 
of [A(Mgb)/A(log(R))]/[A(FUV-y)/A(log(R))] listed in 
Tabled 

In Figure (9] we show the histograms of the Mg h and 
[Z/H] gradients with respect to the UV colours of individ- 
ual galaxies (e.g. A(Mg &)/A(FUV- V)). We present the 
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Figure 8. Same as Figures [6] and [7| but for the stellar velocity dispersion. 
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Figure 9. Histograms of the individual gradients of the UV colour— line index relations for the sample galaxies. The black solid lines 
arc Gaussian fits to all quiescent galaxies. The blue vertical solid lines show the global gradients determined from fitting the luminosity- 
weighted galaxy profiles (solid lines in Figures [l] and |3} , while the dashed lines indicate the associated uncertainties. 



cases of Mgfo and [Z/H] only, showing the tightest corre- 
lations in Figures [6] and [T] These diagrams are useful to 
test whether the individual galaxy gradients (local gradi- 
ents) follow the global correlations (solid lines in Figures [l] 
and [3| . The black histograms represent UV-weak and U V 
upturn galaxies, whereas the red shaded histograms show 
UV upturn galaxies only. We also fit Gaussians to the UV- 
weak and UV upturn galaxy gradient distributions and plot 



them as black solid lines in each panel. The blue vertical 
solid lines show the global gradients (solid lines in Figures [1] 
and |3]) , while the blue dashed lines indicate the associated 
uncertainties. 

If the global correlations (i.e. between galaxies) also ap- 
ply locally within galaxies, the means of the Gaussian fits 
should be consistent with the global gradients. As discussed 
in Section 13.11 our global correlations are consistent with 
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Figure 10. Same as Figure [9] but for tlic UV coIour-stcUar velocity dispersion relations. The blue vertical solid lines show the global 
gradients determined from fitting the luminosity- wreighted galaxy profiles (solid lines in Figure [4]|. 



those of iPaper XVIlj (determined from a fit to i?c/2 aper- 
ture values for quiescent galaxies). For example, the slopes 
of the best-fit (FUV-NUV)-Mg6, (FUV-NUV)-Fe5015 
and (FUV— NUV)— H/3 relations based on integrated prop- 
erties are -0.64±0.12, -0.21±0.17 and 0.20±0.08 A mag"\ 
whereas our global slopes are -0.73±0.04, -0.36±0.06 and - 
0.19±0.02 A mag , respectively. We note that the distribu- 
tions of the individual galaxy gradients are consistent with 
the global gradients within the errors, although the distribu- 
tions are considerably broader (Figure O. For example, the 
global gradient for (FUV-F)-Mg fe is -0.48±0.03 A mag"\ 
whereas the mean of the Gaussian fit to the local gradient 
distribution is —0.46 A mag~^, with an rms of 0.21 A mag 
Comparing the black and red shaded histograms, one can 
see that the distributions of the UV upturn galaxies are 
narrower than those of the UV-weak galaxies. 

Similarly to Figure[9l Figure fTOl shows the histograms of 
the stellar velocity dispersion gradients with respect to the 
UV colours of UV-weak and UV upturn galaxies. In the case 
of the velocity dispersion, however, the global gradients (ver- 
tical solid lines) are different from the local gradients (i.e. 
individual galaxies). This implies that the UV strength is 
not highly dependent on structural or dynamical properties. 



3.4 Multi-dimensional relations 

Figures [TH31 El and [7] show that the (FUV-NUV)-Mg6 re- 
lation is the tightest correlation. We now consider multiple- 
parameter relations analog o us to the fundamen tal plane 
iDiorgovski fc Davisl 1 19871 : iDressler etal] [l983). Three- 
parameter relations are not trivially expected in the absence 
of a priori physical foundations, such as the virial theorem 
in the case of the Fundamental Plane. To minimise biases 
arising from the very different definitions and hence dynamic 
ranges of the parameters considered (e.g. line strengths, stel- 
lar population properties, velocity dispersion and colours), 
we normalize each quantity to yield an observed range of 
1.0. For example, the observed range of FUV— V is roughly 
4.7 to 7.3 mag (a difference of 2.6 magnitudes), and so in the 
abscissa of Figure [TT1 0.0 means 2.6 mag bluer (UV stronger) 
than 1.0. Likewise, for the ordinate of the bottom-left plot, 
the observed range of Mg b is roug hly 3.2 to 5.3 A (a differ- 



ence of 2.1 A), and that of log (tr/kms^^) is roughly 1.9 to 
2.5 (a difference of 0.6). 

We have explored a number of parameter combinations 
and present some cases showing tight correlations in Fig- 
ures [TT] and [13 

The correlation coefficients reported in the bottom-right 
corner of each panel in Figures [11] and [12] are not much bet- 
ter (i.e. larger) than the simple one-parameter fits shown 
in Figures [1] and [3] This means that the relative FUV 
strength is not significantly better described when a third 
param eter is employ ed. This is in agreement with our find- 
ing in jPaper XVIIll that it is stellar population properties 
rather than galaxy dynamical properties that govern the UV 
strength. 



4 THE UV-INDEX RELATIONS 

In summary, the results presented in Section [3] indicate tight 
local and global correlations between UV colours and Mgfc 
and marginal correlations between UV colours and [Z/H]. 
These UV-Mg6 and UV-[Z/H] relations are all the more 
significant because the global gradients are consistent with 
the local gradients. If the UV strength is related to the Mg 
line strength, galaxies with larger Mg gradients should have 
larger UV colour gradients, and the global gradient across 
galaxies should be consistent with the local gradients of in- 
dividual galaxies as observed. Therefore, our results suggest 
that Mg is a main driver of the UV upturn or Mg and the 
UV upturn share the same primary driver, and the correla- 
tions between UV colours and [Z/H] can be understood by 
the role of Mg in the [Z/H] measurement. 

Mg being the most representative a element, it is puz- 
zling that [a/Fe] does not show a cl ear correlation with 
UV. Meanwhile, it should be noted that lCarter et al.l (|201l[ ) 
found good correlations between integrated FUV— NUV 
colours and integrated [a/Fe] and [Z/H]. The weak corre- 
lations between the FUV colours and [a/Fe] in this paper 
may be due to the fact that the a-element is not directly 
measured from line strengths but derived from model fitting 
and so much more difficult than measuring a single metal 
line (e.g. Mg6). Furthermore, the stellar population prop- 
erties ([a/Fe], [Z/H] and age) depend on the stellar models 
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Figure 13. Radial profUcs of ago and mctallicity ([Z/H]) for UV upturn galaxies. The values are derived from NUV-defined elliptical 
annuli at JJc/4, i?c/2 and 2/?e/3 . Each profile is given a slightly different x'-axis offset for clarity. 




Figure 14. Same as Figure [131 but for UV-weak galaxies. 



and spectral librar y adopted (see e.g. ICoelho et al] l2007l : 
iPercival et"aLll2009t ). 

Given the reasonably-clear correlation between [Z/H] 
and UV colours, the apparently-weak correlation shown by 
Fe5015 and UV colours seems mysterious. A possibility is 
that Fe5015 may not be a robust tracer of [Z/H] because 
it is heavily conta minated by poor ly understood absorption 
lines in this band (|Lee et al.ll2009l ). 

We have also shown that global correlations exist be- 
tween UV-related colours and the H/3 line strength, but we 
find no correlation between UV-related colour gradients and 
internal H/? gradients. This can also be explained by mctal- 
licity effects on H/? as discussed in Section 3.1. 

Given the above, the metallicity gradients of UV up- 
turn galaxies must be preserved to t he present day. Is this 
possible in the hierarchical paradigm? iHopkins et al.l (|2009l ) 
carried out simulations and claimed that metallicity gradi- 
ents are only weakly af fected by both wet and dry mergers. 
iDi Matteo et al.l (|2009l ') also argued that metallicity gradi- 
ents can be preserved in major dry mergers if the metallicity 
gradient of companion is sufficiently steep. 



Lastly, we point out that galaxies having recent star 
formation are outliers in the UV-Mgfe and UV-[Z/H] rela- 
tions; they usually have lower Mgfo and [Z/H] values. It is 
therefore clear that tight UV-Mgfe and UV-[Z/H] relations 
are only present for quiescent galaxies. 



5 ORIGIN OF THE UV UPTURN 

5.1 Stellar population parameters from line 
strengths 

If different UV strengths are caused by different stellar pop- 
ulation properties, this implies a difference between the stel- 
lar populations of UV upturn and UV-weak galaxies. To at- 
tempt to understand the conditions that create the observed 
UV upturns, we explore the stellar population parameters 
after separating quiescent galaxies into UV upturn and UV- 
weak galaxies, based on their UV radial colour profiles as 
discussed in Section 3.3. 

Figure [13] shows the age and [Z /H] values of U V up- 
turn galaxies at three different radii. To derive these stellar 
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Figure 15. The moan radial age and metallicity p rofiles of UV upturn (blaclc solid linos) and UV-woak (red dotted lines) galaxies based 
on the derived stellar population parameters from lPaper XVl3 . 



population parameters, we used here NUV-defined elliptical 
annuli of semi-major axes i?c/4 (i.e. central), i?c/2 (i.e. inter- 
mediate) and 2Rc/3 (i.e. outskirts), respectively. For clarity, 
each profile is given a slightly different a;-axis offset. In UV 
upturn galaxies, the derived ages extend from roughly 11 to 
17 Gyr. In the case of [Z/H], the range is roughly —0.25 
to -f 0.05. We note that the age gradients are generally flat, 
whereas significant metallicity gradients exist. 

Similarly to Figure 1131 Figure [14] shows the age and 
[Z/H] gradients of UV-weak galaxies. The age gradients are 
again nearly flat, even though the metallicity gradients are 
weaker than those of UV upturn galaxies. 

For comparison, we construct the mean radial age and 
metallicity profiles for the samples of UV upturn and UV- 
weak galaxies. The results are shown in Figure [151 fhe er- 
rors correspond to the standard deviations about the mean. 
There is a notable difference in age between the two sam- 
ples (left panel), and a slight [Z/H] difference (right panel), 
especially in the central regions. Based on this, one might 
think that it is the age difference that governs the UV 
strength. However, without a radial age gradient, age alone 
cannot explain the central concentration of FUV light ob- 
served. Metallicity gradients are steep, especially in UV 
upturn galaxies, and so are probably essential to explain 
the centrally-concentrated FUV excess (as suggested in Sec- 
tions [3] and [4]) . 



5.2 New stellar population parameters from 
colours 

If it is indeed stellar population properties that govern the 
UV strength of a quiescent galaxy, as we concluded in Sec- 
tions [3] and [?] then the same stellar population models 
should be able to reproduce the UV relations for both in- 
tegrated and internal properties. To find the range of age 
and metallicity that reproduce the observed UV properties, 
we attempt to fit the observed colour profiles (rather than 
line strengths) with the stellar population synthesis mod- 
els of 1x1 (|2003l ) and lYi. Demargue fc Oemlerl (|l997l . Il998l ) 
that were calibrated to reproduce the observed properties 
of the UV upturn phenomenon. These models are sensi- 



tive to metallicity, age, stellar mass loss, helium enrichment 
(AY/AZ) and the distribution of stars on the HB (see e.g . 
Dorman. O'Connell fc Roodlll995l :[ Percival fc Salarisll201ll : 



Chung. Yoon fc Leell20l'li r 



The comparisons between the observed colours and 
models are shown in Figures [16] and 1171 The existence of 
strong radia l gradients of UV colours within galaxi es is well 
known (e.g. iO'Connell et al.|[l993 : lOhl et al.lil998l ). and in 
early-type galaxies they are likely the result of underlying 
stellar population gradients, mainly in age or metallicity or 
both. Therefore, we consider two simple cases: (1) an age 
gradient and (2) a metallicity gradient. 

In the case of the metal-rich model, suggesting 
that the metal-rich HB stars are the most likely 



UV sources (se e Section [T] Greggio fc Renzinil 



Horch et all! 199 2*: 'Brcs sm et al.lll994l : iDorman et all 



199C 



1995 



Yi. Demargue fc Kim 1 9971 ). we could not reproduce the ob- 



served colour gradients at all using an age gradient alone. 
However, if we consider an internal metallicity gradient 
(with no age gradient), metal-rich models are able to repro- 
duce the observed colour profiles. We assume, therefore, that 
the internal colour gradients of quiescent early-type galaxies 
are the result of simple metallicity gradients. The overall age 
of an individual galaxy is permitted to vary over the range 
1 to 20 Gyr. 

To constrain the age and metallicity of each galaxy, we 
thus fit the model to the observed FUV, NUV and V surface 
brightnesses at the same three radii used previously (center 
{Re/4:), intermediate radius {Rc/2) and outskirts (2J?e/3)), 
and compute the associated statistics to obtain a prob- 
ability distribution of the age and metallicity. Figure 1161 
shows the observed FUV— NUV and FUV— V colour profiles 
(filled symbols) and the best-fit model profiles (solid lines) of 
NGC4552 (blue triangle), 4278 (black circle) and 4477 (red 
diamond) . Three open symbols of each galaxy represent the 
three fitted regions. This figure shows that we have found a 
plausible range of metal-rich model parameters that is capa- 
ble of reproducing the observed amplitude of the UV upturn 
with the conventional values of AY/AZ = 2.5, -q — 0.7 and 
[a/Fe] = 0.0. 

Metal-poor models, suggesting that the dominant UV 
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Figure 16. Comparison between the observed colours and metal-rich models of three sample galaxies (NGC4552, 4278 and 4477). The 
FUV— NUV (left) and FUV— V [right) colours of each galaxy are shown as a function of radius . To derive the age and metallicity at a 
given radius, we fit the observed colours with the stellar population synthesis models of lYil l|2003l l. Three different regions for each galaxy 
are used (open symbols; center, intermediate radius and outskirt). The best- fit model is shown as a solid line. 
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Figure 17. Same as Figure [T6l but for the metal-poor model. 



sources are older and me tal-poor HB star s and their post-HB 
progeny (see Section [Tl iLee et all 1 19941 : IPark fc Le3ll997l : 
iBuzzoni fc Gonzalez-Lopezliral 20081 ). fail to reproduce the 
observed surface brightnesses of UV upturn galaxies (e.g. 
NGC 4552 and 4278) when we consider age and metallicity 
gradients, as shown in Figu re II 7[ although th is depends on 
the adopted models (see e.g. lPark fc LedFigOTl ). We note that 
the model-fit for UV weak galaxy (e.g. NGC 4477) based on 
metal-poor models is also reasonably good. A number of 
recent studies found peculiar globular clusters having EHB 
stars that can be explained by the presence of super- helium- 
rich populations (|Lee et al.ll2005. , .20071), which could also be 
a copious source of FUV flux. Indeed, I Chung et al.l (|201ll ) 
claimed that their models with helium-enhanced subpop- 
ulations can reproduce the observed spectra of NGC 4552 
and 4649. To perform the same test for this metal-poor but 
helium enhanced hypothesis, however, require super-helium- 
rich stellar population models which are not fully available 
yet. 



Analogously to Figure [T3l the best-fit values of age and 
metallicity based on metal-rich models, for all galaxies clas- 
sified as UV upturn from their radial colour profiles, are 
shown in Figure (TH] For clarity, the profiles are again given 
slight offsets with respect to each other along the s-axis. Er- 
rors show parameter values within Icr of the best-fit model. 
According to these models, the age range is 9 to 12 Gyr and 
the metallicity range is roughly —0.05 to +0.28. Similarly 
to Figures [T4l and [181 Figure [191 shows best-fit values of age 
and metallicity for UV-weak galaxies. In this case, the range 
of ages is 9 to 11 Gyr, with metallicities of roughly —0.27 to 
+0.20. 

Like Figure [151 Figure [20l shows the mean stellar pop- 
ulation properties and the standard deviations about the 
means for UV upturn and UV-weak galaxies (filled squares) . 
According to the models, the mean stellar age of UV-upturn 
galaxies (black solid lines) is about 11 Gyr and their cen- 
tral regions have a metallicity distribution peaking around 
0.2, while at the effective radius the metallicities decrease to 
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Figure 18. Radial profiles of age and m etallicity for UV upturn galaxies, as derived from fits to the observed colours and the stellar 
population synthesis models of|Y| ||2003| ) at -Rc/4, -Rc/2 and 2Rc/3. Each profile is given a slightly different x-axis offset for clarity. 
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Figure 19. Same as Figure [TSl but for UV-weak galaxies. 
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Figure 20. The mean radial age and met allicity profiles of UV upturn (black solid lines) and UV-weak (red dotted lines) galaxies based 
on colour fits and the models of IyH liooi) (filled squares). For comparison, we also show the mean ages and mctallicities of UV upturn 
and UV-weak galaxies based on the derived stellar population parameters of lPaper XVIll (filled circles; see Figure [TSll . 
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slightly sub-solar. UV-weak galaxies (red dotted lines) have 
a similar mean stellar age but significantly smaller metallic- 
ity gradients compared to UV upturn galaxies. The metal- 
licity difference between UV upturn and UV-weak galaxies 
in the central regions produces the difference in the UV up- 
turn strength. The UV upturn phenomenon seems indeed to 
be due to stellar population properties. 

For comparison, we plot again the mean ages and metal- 
licities derived from spectral line strengths (see Figure fTS)) as 
filled circles. In both cases, UV upturn galaxies are slightly 
older than UV-weak galaxies and the metallicities of UV up- 
turn and UV-weak galaxies reveal differences, especially in 
the central regions. The trends of age and metallicity derived 
from spectral line strengths and colours are thus similar, but 
the actual values reveal significant differences. This can be 
explained in terms of the age-meta llicity degene racy. 

The previous SAURDN work of IPaper XVlj found that 
models with very high ages (»il4 Gyr) match the high 
Mg absorption of UV upturn galaxies. In contrast, through 
our colour fitting method, we find that slightly younger 
(«11 Gyr) but more metal-rich models better reproduce the 
pronounced UV fiux in the interiors of UV upturn galaxies. 
It is well known that line strengths and colour indices suffer 
from the age-metallicity degeneracy (Worthey 1994). Fur- 
thermore, as described in Section [2.31 we should not believe 
the absolute ages and metallicities returned by the models, 
even though they still allow us to make relative comparisons 
between gal axies. The m odel parameters derived from line 
strengths in lPaper XVlj (Section lS.ip and from colours here 
(Section 15. 2p can thus probably be reconciled through the 
age-metallicity degeneracy. 

Therefore, despite having failed to find the same param- 
eter ranges when using spectral line strengths and colours, 
we do feel that we have found plausible parameters that re- 
produce the observed radial colour profiles of UV upturn 
and UV-weak galaxies. Furthermore, it is important to re- 
produce the UV relations for internal pr o pertie s using the 
stellar population synthesis models of|Yj (|2003h . which are 
specifically designed to explain the UV upturn phenomenon. 
Notably, our models reproduce the range of UV strengths 
among elliptical g alaxies by varying the metallicity. In the 
context of the |y| (|2003l ) models, the UV upturn is a stel- 
lar phenomenon associated with the metal-rich (high Mg6 
line strength) population of early-type galaxies, and the UV- 
bright central regions of giant elliptical galaxies (e.g. NGC 
4552 and 4278) are enhanced in metals by roughly 60% com- 
pared to the UV-weak parts. 



6 SUMMARY 

We have used ultraviolet photometry from GALEX, ground- 
based optical photometry from MDM and ground-based op- 
tical integral-field spectroscopy from SAURQN to study the 
UV-line strength and UV-stellar population relations of 
34 early-type galaxies from the SAURON sample. Our main 
advantage in exploring these correlations is that all photo- 
metric and spectroscopic measurement are derived using the 
same ellipses and apertures. 

First, we have found that passively-evolving galaxies 
show clear correlations between the colours involving FUV 
{FUY-V and FUV-NUV) and the Mg b and H/3 absorption 



line strength indices, as well as the metallicity [Z/H], and 
that these correlations are strongest for FUV-NUV. This 
result is supported by a pixel-by-pixel analysis. This suggests 
that either Mg is a main driver of the UV upturn or Mg and 
the UV upturn share the same primary driver, and that 
the correlation between the UV colour and metallicity is 
explained by the role of Mg in the metallicity measurement. 
Likewise, the correlation of the UV colour with H/3 can be 
explained by the metallicity effect on H/3. 

In addition, we have also derived logarithmic internal 
colour, measured line strength and stellar population gradi- 
ents for each galaxy and again found a strong dependence 
of the FUV-V and FUV-NUV colours on the Mgfe line 
strength, as well as on metallicity. Moreover, the Mg6 and 
[Z/H] global gradients with respect to the UV colours (e.g. 
A(Mg& )/A(FUV— V)) are consistent with their respective 
local gradients, implying that the global correlations also 
hold locally within galaxies. 

Finally, if the UV strength is indeed governed by stellar 
population properties, then the same models should be able 
to reproduce the UV relations for both integrated proper- 
ties and internal properties (within each galaxy). From sim- 
ple models based on UV— optical colour fits of UV upturn 
and UV-weak galaxies, we have identified a plausible range 
of parameters that reproduces the observed radial colour 
profiles of each galaxy type. In these models, the UV up- 
turn is a stellar phenomenon associated with the metal-rich 
(high Mg6 line strength) population of early-type galaxies, 
and the inner parts (UV-bright regions) of giant elliptical 
galaxies are enhanced in metals by roughly 60% (e.g. NGC 
4552 and 4278) compared to the outer parts (UV-weak re- 
gions). We stress here that stars originating the UV upturn 
phenomenon are detected in the central region of galaxies 
(especially within Re/4), and therefore they are contribut- 
ing < 20 % of the total stellar light. This means that the UV 
upturn phenomenon is naturally produced if the galaxy has 
a small fraction of Mg enhanced populations. Assuming that 
the galaxy has a 10% metal-enhanced population, it cannot 
make a significant change in the overall metallicity of the 
galaxy, but this population can evolve to the FUV excess 
stars. 

Without doubt, we now have much more information 
on the UV upturn than a decade ago. However, there still 
are many puzzles. For example, it is surprising that the a- 
enhancement ([a/Fe]) does not show a correlation with UV 
strength, considering that Mgb shows a tight correlation 
with UV colours and [Z/H] also shows a marginal corre- 
lation. Improvements over the work presented in this paper 
will come from a more careful population modelling analysis. 

Our achievements here were mainly obtained through 
the spatially-resolved analyses that have just recently be- 
come possible. A critical weakness of our work however is 
the poor spatial resolution of the UV data. Breakthroughs 
will be guaranteed when higher spatial resolution (arcsec 
scale) FUV data become available, hopefully in the near fu- 
ture. 
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